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ABSTRACT 

Continuous and precise space-based photometry has made it possible to measure the orbital 
frequency modulation of pulsating stars in binary systems with extremely high precision over 
long time spans. Frequency modulation caused by binary orbital motion manifests itself as a 
multiplet with equal spacing of the orbital frequency in the Fourier transform. The amplitudes 
and phases of the peaks in these multiplets reflect the orbital properties, hence the orbital 
parameters can be extracted by analysing such precise photometric data alone. We derive an¬ 
alytically the theoretical relations between the multiplet properties and the orbital parameters, 
and present a method for determining these parameters, including the eccentricity and the ar¬ 
gument of periapsis, from a quintuplet or a higher order multiplet. This is achievable with the 
photometry alone, without spectroscopic radial velocity measurements. We apply this method 
to Kepler mission data of KIC 8264492, KIC 9651065, and KIC 10990452, each of which is 
shown to have an eccentricity exceeding 0.5. Radial velocity curves are also derived from the 
Kepler photometric data. We demonstrate that the results are in good agreement with those 
obtained by another technique based on the analysis of the pulsation phases. 

Key words: asteroseismology - stars: oscillations - stars: variables - stars: binaries - stars: 
individual (KIC 8264492; KIC 9651065; KIC 10990452) - techniques: radial velocities. 


1 INTRODUCTION 

The study of binary and multiple stars provides the best mea¬ 
surements of stellar parameters, giving direct, nearly model- 
independent knowledge of stellar masses, and radii as well from 
eclipsing binaries. These are fundamental to our understanding of 
stellar structure and evolution, the bedrock upon which much of 
stellar astrophysics rests. Asteroseismology can now also provide 
stellar masses and radii, and even ages for individual stars, but our 
confidence in asteroseismology rests on calibration of asteroseis- 
mic masses and radii with those determined from binary stars. 

The primary observational data for binary star studies are the 
light curve and radial velocity curve. Additional constraints are 
given by measurement of stellar effective temperature, T e g, from 
spectroscopy, and, for a few stars, stellar radius from interferome¬ 
try. In the case of T e s and interferometric radius, single measure¬ 
ments are sufficient. The light curve and radial velocity curve, how¬ 
ever, are by definition time series, and it is these that demand so 
much time and effort from observers. 

Traditionally, light curves have been measured with ground- 
based photometry from a single observing site in one or more pho¬ 


tometric bandpasses. Radial velocity curves have been measured 
by obtaining spectra from a single site repeatedly until the orbital 
phase is well covered. So, a large amount of observing time is re¬ 
quired, often on large telescopes for fainter stars. This is the tradi¬ 
tional bottleneck to binary star studies, and until recently the result 
has been that high precision masses and radii for stars were only 
known for a matter of hundreds of stars. That all of stellar structure 
and evolution theory, and beyond that galactic astronomy, should 
be based on such a sparse base has been a significant concern. 


Now the study of stellar light curves has been revolutionised. 
Whereas until recently light curves were studied for no more than 
a few stars from a few sites, for a limited number of nights per year 
with great gaps in the data, we now have nearly continuous space- 
based photometry for hundreds of thousands of stars. Whereas from 
the ground photometric measurements were precise to parts per 
thousand, or in the best cases, about 10 times better than that, the 
space data give precision as good as parts per million. The satellite 
telescope photometry revolution has a pedigree from the Canadian 
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MOST missioi 
Kepler missiorU 


through the ESA CoRoT missior@ to the NASA 

_ and it anticipates the upcoming NASA mission 

TESS0and ESA mission PLATC0. 

It is the Kepler mission that has led to the greatest results, 
with - as of this writing - over 4600 exoplanet candidate discov- 
eries^llHfl over 1800 confirmed exoplanets, asteroseismic results 
for thousands of pulsating stars, and with light curves of more 
than 2700 eclipsing binary starfl The main Kepler mission data 
set provides 4 yr of nearly continuous photometric white light data 
with 30-min time resolution for 150 000 stars, with further data for 
shorter time spans for another 50 000 stars. The follow-up studies 
for exoplanets involve many hundreds of astronomers and require 
large amounts of ground-based telescope time - often on some of 
the world’s largest telescopes - for radial velocity measurements 
necessary to confirm that masses are indeed planetary. The thou¬ 
sands of binary stars also require large amounts of ground-based 
telescope time to acquire the radial velocities needed for full astro- 
physical study of the multiple systems. 

The impediment in these studies is the acquisition of the spec¬ 
tra that have traditionally been used to measure radial velocity from 
the Doppler shift of wavelength of spectral lines. However, Doppler 
shift can be measured for any astrophysical source of a stable fre¬ 
quency. We have developed two methods to do this directly from 
Kepler mission photometric light curves for pulsating stars, with¬ 
out recourse to ground-based observations, where the pulsation 
frequencies are the standard, rather than spectroscopic wavelength. 
These methods obviate the need for any ground-based observations 
to determine radial velocity, and they provide nearly continuous 
radial velocity curves that also cover the entire Kepler 4-yr obser¬ 
vational time span - an unreachable goal with ground-based radial 
velocity measurements. Our methods are limited to stars with sta¬ 
ble pulsation frequencies, but there are thousands of these in the 
Kepler data with amplitudes high enough to apply our techniques. 

Our two methods are t he frequency modulation (FM) method 
IShibahashi&Kurtz| 1 20121) and the pha s e modulati on method 
(Murp hy et al. 12014 : see also lBalonall2014l : lKoenll2014l) . Our tech¬ 
niques are related to older observed-minus-calculated (O — C) 
methods of studying stellar frequency variability, but the FM and 
PM techniques both have the advantage that they use the entire 
data set for maximum signal-to-noise ratio and that they are par¬ 
ticularly suited to multi-mode pulsators. The FM technique has the 
highest possible frequency resolution, and the PM technique is both 
easily automated and can combine the results from many pulsation 
frequencies easily. In addition, many studies of pulsating stars for 
asteroseismic inference begin with the frequency spectrum, and the 
FM technique shows explicitly the patterns expected in an ampli¬ 
tude spectrum, or power spectrum, for pulsation frequencies under¬ 
going periodic Doppler shifts from binary orbital motion. 

We have demonstrated the validity of the FM method by show¬ 
ing the consistent results obtained from it when compared t o a tra¬ 
ditional eclipsing binary light curve analysis l iKurtz et al]l2015l) . 


1 http://most.astro.ubc.ca/ 

2 http://sci.esa.int/corot/ 

3 http://kepler.nasa.gov 

4 http://space.mit.edu/TESS/ 

5 http://sci.esa.int/plato/ 

6 http://exoplanets.org/ 

7 http://exoplanetarchive.ipac.caltech.edu 

8 http://kepler.nasa.gov/Mission/discoveries/ 

9 http://keplerebs.villanova.edu 


Furthermore, the FM technique is also applicable to non-eclipsing 
systems and is a powerful tool to determine the orbital parameters. 

In practice, the information desired from the radial ve¬ 
locity curve includes the orbital period, the mass function, 
/(mi, m 2 , sin i), the orbital semi-major axes, ai sin i and a 2 sin i, 
the eccentricity, e, the angle between the nodal point and the peri- 
apsis, zd, and the time of periapsis passage, t p . Shiba hashi & Kurta 
l l2012l) showed how to derive all of these exc ept e and w fro m 
photometric data using the FM method, while lMurphv et IDGoU 
showed how to derive all of them from the PM method. In this pa¬ 
per we show how to extract e and w with the FM method. We also 
show how to derive a radial velocity curve from the FM method, 
because many investigators have developed software to analyse bi¬ 
nary star orbits using light curves combined with radial velocity 
curves. For example, the pre-eminent program for binary star anal¬ 
ysis is now PHOEB^EJ which is designed for radial velocity curve 
input JPrsa & Zwittej2005t) . We emphasise that the FM technique 
can provide directly all the information that is extracted from the 
radial velocity curve, but we also provide the method to generate 
the radial velocity curve itself, both for input to programs such as 
PHOEBE, and because the radial velocity curve is visually helpful 
for thinking about the binary star orbit. 

We have discussed Doppler shift to help visualise what our 
techniques do, but the FM and PM methods do not rely directly on 
Doppler shift; they are built on the equivalent Rpmer time delay. 
In Sections 2-4 of this paper we derive the relations to determine 
the mass function, orbital semi-major axis for a pulsating compo¬ 
nent, eccentricity, w and t p , as well as the radial velocity curve, 
from the frequencies, amplitudes and phases of the components of 
a multiplet in the Fourier spectrum, which are split by the orbital 
frequency of the binary star. 

In Section 5, we provide some examples using Kepler data 
for stars with high eccentricity. We compare the results with those 
obtained with the PM method, and discuss the circumstances for 
which FM or PM is the preferable technique. We show that the 
method is robust, and that it is capable of detecting companions 
with masses down to brown dwarf (m ^ 0.08 Mg) and even super- 
Jupiter masses (m SC 0.03 Mg). 


2 THE LIGHT TRAVEL TIME EFFECT ON PULSATION 
IN A BINARY STAR 

Let us consider a star sinusoidally pulsating with a single angu¬ 
lar frequency wo in a binary. We name the stars T and ‘2’, and 
suppose that star 1 is pulsating. As the pulsating star moves in its 
binary orbit, the path length of the light between us and the star 
varies, leading to a periodic variation in the arrival time of the sig¬ 
nal at Earth. The difference in the light arrival time, r(t), compared 
to the case of a signal arriving from the barycentre of the binary 
system is given by 

r(f) := - / itrad (t')dt', (1) 

c Jo 

where c is the speed of the light and « ra d(f) denotes the radial ve¬ 
locity, due to the orbital motion, of the star 1 at the time t, where 
the epoch is the time at which the star passes the nodal point (‘N’ in 
Fig.QJ directed toward us. We follow the convention that the sign 
of the radial velocity is positive when the star is receding from us. 
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Hereafter we call r(i) the “time delay”. With the help of equation 
0 for the time delay, the observed luminosity variation at time t is 
written as 

AL(t) = A cos {ix>o[t — r(t) ]+</} , (2) 

where A is the amplitude of pulsational luminosity variation and f 
denotes the pulsation phase at t = 0. In the following we derive 
the Fourier transform of the luminosity variation with the form of 
equation 0. 


2.1 The radial velocity as a function of time 


First, we have to derive the radial velocity as a function of time. 
Let us suppose a plane that is tangential to the celestial sphere on 
which the barycentre of the binary is located, and let the z-axis that 
is perpendicular to this plane and passing through the barycentre 
of the binary be along the line-of-sight toward us (see Fig0. The 
orbital plane of the binary motion is assumed to be inclined to the 
celestial sphere by the angle i, which ranges from 0 to n. The orbital 
motion of the star is in the direction of increasing position angle, 
if 0 ^ i < 7t/2, and in the direction of decreasing position angle, 
if 7t/2 < i ^ 7t. We write the semi-major axis and the eccentricity 
of the orbit as ai and e, respectively. Let w be the angle measured 
from the nodal point N, where the motion of the star is directed 
toward us, to the periapsis in the direction of the orbital motion of 
the star. Also let / be the ‘true anomaly’, which is the instantaneous 
angle measured from the periapsis to the star in the direction of the 
orbital motion of the star , and let r be the distance between the 
barycentre and the staFl 

To make some complicated formulae derived later in this pa¬ 
per more easily traceable, we must repeat the fundamental deriva¬ 
tion of the radia l velocity as a funct i on of time, although it was 
already given in Ishibahashi & Kurtzl J2012I) . The z-coordinate of 
the position of the star is written as 

z = r sin(/ + vj) sini, (3) 


and the radial velocity, v ra d := — dz/df, is 


^rad — 


dr d / 

sin(/ + vj) + 1 —jj- cos(/ + vj) 


dt 


(4) 


With the help of the know n laws of motion in an ellipse (e.g., 
iBrouwer & Clemencelll96ll) . 

df aifi(l + ecos /) 

r a - Vl -e> (5) 


and 

dr ai fie sin/ 

dt = Tf3?’ (6) 

the radial velocity of star 1 along the line-of-sight is expressed as 

tlai sin i r , r \ n 

v ra d = - , [cos(/ + U7) + e cos vo \, (7) 

V1 - e 2 

where Q denotes the orbital angular frequency. The time depen¬ 
dence of radial velocity is implicitly expressed by the true anomaly 


11 The range of i, the definition of vj, and the sign of z are often differently 
given by different authors. It should also be noted here that we cannot dis¬ 
tinguish a motion illustrated in Fig.|7]and a motion in the inverted image of 
the bottom panel of Fig.[7]from the radial velocity measurement alone. The 
former is a prograde motion, while the latter is a retrograde motion with the 
same value of vj, but with the inclination angle being the supplementary 
angle of i shown in the top panel of Fig.|T| 



Figure 1. Top: Schematic side view of the orbital plane, seen from a 
faraway point along the intersection of the orbital plane and the celestial 
sphere, NFN', where the points N and N' are the nodal points, respectively, 
and the point F is the barycentre of the binary system; that is, a focus of the 
orbital ellipses. The orbital plane is inclined to the celestial sphere by the 
angle i, which ranges from 0 to n. In the case of 0 ^ i < 7t/2, the orbital 
motion is in the direction of increasing position angle of the star, while in 
the case of n/2 < i ^ 7t, the motion is the opposite. The z-axis is the 
line-of-sight toward us, and z = 0 is the plane tangential to the celestial 
sphere. Bottom: Schematic top view of the orbital plane along the normal 
to that plane. The periapsis of the elliptical orbit is P. The angle measured 
from the nodal point N, where the motion of the star is directed toward us, 
to the periapsis in the direction of the orbital motion of the star is denoted 
as vj. The star is located, at this moment, at S on the orbital ellipse, for 
which the focus is F. The semi-major axis is a\ and the eccentricity is e. 
Then OF is aie. The distance between the focus, F, and the star, S, is r. 
The angle PFS is ‘the true anomaly’, /, measured from the periapsis to the 
star at the moment in the direction of the orbital motion of the star. ‘The ec¬ 
centric anomaly’, u, also measured in the direction of the orbital motion of 
the star , is defined through the circumscribed circle that is concentric with 
the orbital ellipse. 


/, which can be written in terms of ‘the eccentric anomaly’, u (see 
Fig 0 , defined through the circumscribed circle that is concentric 
with the orbital ellipse, satisfying 

r cos f = ai (cos u — e) (8) 

and 

r sin / = ai y/l — e 2 sin u. (9) 

The eccentric anomaly u is written as 

u = Q(t — i p ) + 2 —J„(jie) sin nfi(f — f p ), (10) 

n= 1 

where J n (x) denotes the Bessel function of the first kind of integer 
order n, and f p denotes the time of the periapsis passage of the star. 

The trigonometric functions of the true anomaly / are ex¬ 
pressed in terms of a series expansion of trigonometric functions of 
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Figure 2. e(l-e 2 )- 1 / 2 J4 (ne) / Jn (ne) for n = 1 (red), n = 2 (green), 
n = 3 (blue), and n = 4 (magenta) as functions of eccentricity e. This 
shows that e(l — e 2 ) -1 / 2 J' n (ne )/ J n (ne) for n = 1, ...,4is ~ 1 except 
for cases of large e, which means i9(e, to - ) ~ zu except for cases of large e. 
This characteristic is used later in Sect. EH 


Figure 3. The coefficients of £i (e, zu) (purple), £2 (e, ro) (green), £3 (e, zu) 
(blue), and £ 4 (e,zu) (magenta) for the Fourier amplitudes of the phase 
modulation as functions of eccentricity e. The band of each curve shows 
the range of zu from 0 to 2n. The tu-dependence of £ 2 * £ 3 , and £4 is sub¬ 
stantially weaker than £1 in a wide range of e. This figure is used later in 
Sect.l331 


the time after the star passed the periapsis with Bessel coefficients 
(see Appendix 1): 

cos f = — e+— --—V J n (ne) cos nQ,{t — f p ), (11) 

e f 

n =1 
oo 

sin f = 2\/T~~e 2 ^^ J n '(ne) sinnfl(t. — t p ), (12) 

n=1 

where J„'(x) := dj n (x)/dx. The radial velocity is then written 
explicitly as a function of time: 


OO 

ttrad(i) = —flat sin*y n£ n (e,iu) cos [nfi(i — t p ) + #„], (13) 

n= 1 


where 


J 1 _ e 2 

£ n (e,zu) := 2—- J n (ne) 

ne 


x \ I ms 2 777 4- ( — — C I' ne | ) sin 2 w 

\s/T^Jn{ne)) 


and 

i ?n(e, w) := arctan 


e J’ n {ne) 
x/1 — e 2 Jn(ne) 


tantr7 


(14) 


(15) 


where arctan(z) returns the principal value of the inverse tangent 
of x. 

Since 


J' n {ne) _ 1 ne/2 n 2 e 2 /2 n 3 e 3 /2 

J n (ne) e (n + 1) - (n + 2) - (n + 3) - " ' 


(16) 


where the series in square brackets denotes continued fractions, the 
^-dependence of £„(e, ur) is weak for e -C 1, but it becomes 
substantial with the increase of e (see Fig.[3]l. Also, since 


MX) = y (-1)* 

k =0 


(x/2) n+2k 
(n + k)\ k\ ’ 


(17) 


the lowest order of the series expansion of (e, zu) in terms of 
power of e is e n ~ 1 ; £„(e, zu) ~ 0(e"~ 1 ) for e«l. 


2.2 The time delay as a function of time 

Integrating equation J 1 31) . we obtain an expression for the time de¬ 
lay as a function of time 

. . OO 

r(t ) = - 0,1 sm 1 y Cn sin[nf2(f - f p ) + #„] - to, (18) 

C n—1 

where 

. . OO 

, v sin i \—> . . , N 

ro(e, zu) :=-\ £„ sm(—nfif p + #„). (19) 

c ' 

7i=l 

Fig.[3]shows £i, £ 2 , £ 3 , and ^(e, ro) as functions of eccentric¬ 
ity e. The band of each curve shows the range of vj from 0 to 27t. 
It is seen that the ro-dependence of £ 2 , £ 3 , and £4 is substantially 
weaker in a wide range of e than is the case of . This is because 
the terms £ 2 , £ 3 , £4 and the higher terms represent a departure from 
a pure sinusoidal variation with the orbital angular frequency Q, 
which is essentially caused by the eccentricity of the orbit. On the 
other hand, £i(e,ro) shows substantial ro-dependence, in partic¬ 
ular with increase in e. This is because the projected light travel 
time is dependent not only on the eccentricity of the orbit, but also 
on the direction of the apparent elongation of the projected orbit, 
particularly in cases of high eccentricity. 


2.3 Phase modulation 

Ishibahashi & Kurtzl I 2 OI 2 I 1 define a parameter a that measures the 
amplitude of the phase modulation when the pulsation frequency is 
treated as fixed: 


a := 


aicog sin i 


( 20 ) 


which is the ratio between the light travel time across the projected 
semi-major axis and the pulsation period of the mode in consider¬ 
ation. It should be noted that the value of a is dependent on the 
mode frequency. The larger the orbit size and the shorter the pul¬ 
sation period, the larger the resulting phase modulation. Then the 
observed luminosity variation, whose amplitude is assumed to be 
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unity, is rewritten as 

OO 

(wot + <p) + a ^2 £n sin (nQt + On) , 

71=1 

where 


A L(t) = cos 


( 21 ) 


<p :=(fr + woto (22) 

and 

On ■= On - ntttp. (23) 

If we regard equation ED as luminosity variation with the intrin¬ 
sic angular frequency cuo, the pulsation phase is regarded as being 
modulated by the orbital motion. 


3 METHODOLOGY FOR DETERMINING BINARY 
PARAMETERS FROM PHOTOMETRY 


The problem we address here is how to reproduce the radial veloc¬ 
ity v Ta d(t) from the observed luminosity variation given by equa¬ 
tion CD. where v Ta d (t) is rewritten as 


^rad(^) — 


Qc 

- a 

coo 


E 


n£ n cos[nflf + 6 n ], 


(24) 


or equivalently, how to deduce the orbital parameters, that is, the 
orbital angular frequency f l, the semi-major axis ai, the eccentric¬ 
ity e, the angle between the periapsis and the nodal point zu, and 
the time of periastron passage of the star f p . 


3.1 Fourier transform of the luminosity variation 

Here we describe the Fourier analysis of the luminosity variation, 
given by equation CD- showing periodic phase modulation. This is 
done with the help of Bessel functions, written with arbitrary real 
x and 9: 


OO 

e ±ix sin 0 _ ^ J n (x)e ±in8 

n= — oo 


(25) 


Applying this, we rewrite equation CD as 


A L(t) = 




x exp 


i(cjQt+<^) 1 T ia:£ n sin(nQ,t-\-9n) 

II e 

71= 1 

OO oo 

iim v' • • • 

—>-oo ‘ ^ ^ 

k\ = — oo fc^r = — oo 

N 

i k n (nQt + 0 n 


(a, 0 t + ¥>) 


Jk n (ot^n) 


(26) 


where 5R(a:) means the real part of x and N denotes a large number. 
This means that a frequency multiplet appears around the intrinsic 
frequency u>o in the frequency spectrum, and that each component 
of the multiplet is separated from its neighbouring components by 
the orbital frequency fl: 


A L(t) = 5R< A 0 e iUot 


+ (-4 +m e i( “ 0+mn)t + 

771=1 


, (27) 


where _4o and A±m are the complex amplitudes of the central peak 
and the sidelobes separated from the central component by ±mfl, 
respectively. 


3.2 Truncation in the case of a < 1 


It is straightforward to compute numerically the terms in equa¬ 
tion 126b for a given set of binary parameters. The purpose of the 
present paper is, however, the reverse: we deduce the binary param¬ 
eters from the Fourier transform of the light curve. An important 
point is that the complex amplitudes, _4o and A± m , are determined 
by the quantities that characterize the binary system, a, £„(e, zu), 
0 n (e, zu), allowing us to inversely determine those binary parame¬ 
ters from the amplitude spectrum of the oscillations. 

In the following, we derive analytical expressions of the 
Fourier transform of the light curve in terms of the binary parame¬ 
ters. In practice, we truncate the infinite expansion series in equa¬ 
tion ( 126b with a finite number of terms, N. There appear (2^+1)^ 
cross-terms, and it is neither practical nor necessary to take account 
of all the terms. Rather, it is better to keep only the leading terms. 
Indeed, in the case of a < 1, the Bessel functions with the argu¬ 
ment of a£ n become negligibly small with the increase of n and 
the order k n . Hereafter, we limit ourselves to consider the cases of 
a < 1; that is, the cases for which the light travel time across the 
projected orbit is longer than the pulsation period. 

Later, in Section 7, we will demonstrate some observational 
examples that show up to septuplets, nonuplets, and even undecu- 
plets. So, here we truncate the infinite series of equation ( 126b with 
N = 5. Then, equation ( 126b produces (2N + 1)^ = ll 5 = 
161051 cross-terms. Among them, we discard terms separated 
from the central peak at w o further than ±6H and keep only the 
undecuplet. For each frequency component of mil, we keep only 
the terms up to 0(a 2 ). 

Then, equation ( 127b shows that the light-time effect in a binary 
star leads to a frequency undecuplet, w = wo ± m H (m = 1,..., 5), 
for which complex amplitudes of the first and second sidelobes are 
given as 


-4±i _, /Ji(a£i)\ iiej 

Ao - UoKi)7 

X f T ( e ±i((9 2 -2fli) 


(28) 


and 


A±2 

Aq 


( Ji(a&) \ 

\Jo(a&)J 


e ±ie 2 


x 


f Ji(aj2) 

\Ja(at,2) 


-1 


( J2(a£i) \ 

\Jo(a(,i)) 


^1(02-20!) 


,(29) 


respectively. The complex amplitudes of the higher order sidelobes 
are given in Appendix 2. 


3.3 Amplitude ratios and phase differences between the 
sidelobes and the central component 

Let Aq and (fro be the amplitude and the phase of the central com¬ 
ponent of the multiplet at the angular frequency wq, and let 
and (fr±m be the amplitudes and the phases of the ±m-th sidelobes, 
respectively. In the following, we derive the amplitude ratios and 
the phase differences between the sidelobes and the central compo¬ 
nent. 
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3.3.1 The case of the first side lobes 


From equation 428k the amplitude ratio between the first sidelobes 
(m = 1) and the central peak is given, up to the order of 0 (a£ i), 
by 


A+i + A -1 n i 

OiKi)' 

Ao 

V Jo(a£i), 

and 


A+i — A i 

( Ji(«6) 

A+i + A_i “ 



cos(20i 


0a). 


(30) 


(31) 


Equation 43 11 1 indicates whether A +i or A -1 has the higher ampli¬ 
tude, according to the sign of cos(20i — 1 ) 2 ). If e <C 1, 0i ~ ot 
( see Sect. l2.lt . Hence, in such a case, whether A +1 or A -1 has the 
higher amplitude is determined by cosot. For 0 ^ ot < n/2 or 
3n/2 < ot ^ 2n, A -1 > A+ 1 , while, for n/2 < ot < 2>n/2, 
A+i > A- 1 . When ot = n/2 or ot = 3n/2, that is, when the 
periapsis or the apoapsis is towards the observer, the amplitudes of 
A +1 and A-\ are equal to each other. 

Also, the phase differences between the first sidelobes and the 
central component are given as 


0±i - 0o := arg ± ^ 

/ 7T 7T 


Let us choose the zero point for the phases so that the phases of 
the first sidelobes are equal. This condition is fulfilled twice during 
one orbital period, at t = to being 


i (fit + # 1 ) + 


Tito + #1 



(33) 


where 


k = 


0 and 1 if 0 < 81 ^ n /2 

1 and 2 if tt/2 < 61 < 3rt/2 (34) 

2 and 3 if 37t/2 ^ 9i < 2n. 


At each of these phases, the Fourier component proportional to 
cos f It of the radial velocity becomes zero. In the case of a circular 
orbit, this corresponds to the moment at which the orbital motion 
of the star is perpendicular to the line-of-sight. At f = to, 

4 >±i{to) — 4>o = ± ± n 

+ (xM) sin(2<,1 -* ) ' <35) 

which means that the phases of the first sidelobes apparently dif¬ 
fer from that of the central peak approximately by either +7t/2 or 
—n/2. This can be used to confirm that FM is caused by the or¬ 
bital motion. Equation 435b also indicates that the deviation of the 
phase difference from +7t/2 or —n/2 is dependent on the sign of 
sin(202 — 0i). Since 0i ~ ot, this means that the phase difference 
is slightly larger (smaller) if the periapsis is in the far (near) side of 
the orbit, with respect to us. 


3.3.2 The case of the second sidelobes 
Similarly, up to the order of 0 (a^ 2 ), 

A + 2 + A -2 ^ 2 f 

Ao \ Jo ( 0 ^ 2 ) / 


( 36 ) 


and 


A +2 — A -2 
A +2 + A_2 


f Ji( q 6) 
V Jo(«^2) 


-1 


( Js(a£i)\ 

u«i); 


cos(20i 


02). (37) 


Equation 437l l indicates whether A +2 or A -2 has the higher ampli¬ 
tude is determined by the sign of 008(2171 — 0 2 )■ Note that if A +1 
is higher (lower) than A- 1 , then A +2 is lower (higher) than A- 2 . 
Also the phase differences are 


4 >± 2 


(j >0 := arg ± 2 Q.t 


~ ±(2m + 6, 2 ) 

( Jija.& y 1 ( J 2 {atji) \ 

\Jo(a&)J \Jo(a£i)J 


sin(20i 


02). 


Hence, at t = to. 


(38) 


0+2(60) — 0—2 (60) — — 7t + 2(2fIfo + #2). (39) 

Combining equations 433b and 439k we obtain 
2(6*2 - 20i) = [0+2(to) - <(>-2(to)] - (4 k + l)7t, (40) 

equivalently, 

201 — 02 = — - [0+2(to) — 0—2(to)] + — (41) 

It should be noted that there remains an uncertainty of 2n7t in the 
measured phase difference, 0 + 2 (to) — 0 - 2 (to), where n is an in¬ 
teger, and then we cannot uniquely determine (20i — 02 ) from 
this equation alone. This uncertainty can be eliminated, however, 
by taking account of the inequality of A± 1 and the deviation of 
0±i(to) — 0o from ±7t/2. 


3.3.3 General description: the higher-order sidelobes 
Similarly, up to 0{ at/m), for the m-th sidelobes (m ^ 1), 

A+m + A — m ^ f 

Ao - VJo(aU); 

Since the value of the left-hand side is observationally obtained, 
this equation determines the value of a£ m . A practical numerical 
method of solving equation 442b is given in Section ldTl The graphic 
solutions are shown in Fig.[4] Also, 

0+m(to) - 0-m(fo) = -7t + 2(mfit 0 + 9 m ). (43) 

This relation allows us to deduce 

2 (9m - md i) = [0 +m (f o ) - 0_ m (f o )] - (m - l)7t. (44) 


3.4 Formal reconstruction of radial velocity from 
photometry 

Equation 442b indicates that the value of a£ m is observationally de¬ 
termined from the amplitude ratio between the sum of the sidelobes 
and the central peak. The graphical solution is shown in Fig.[4] 
The solution is numerically obtained with sufficient accuracy by 
expanding the Bessel functions J n (x) up to the fifth power of x. 
On the other hand, the phase differences between the sidelobes de¬ 
termine 9m from equation 444b . where 9\ is given by equation 433b . 

By substituting the values of a£ m and 9 m into equation 424b . 
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X 



Eccentricity 


Figure 4. The ratio of 2Ji(x)/Jo(x) as a function of x. For x 1, it is 
well approximated by x as shown. For the whole range of 0 ^ x ^ 1, it is 
well approximated by a polynomial of the fifth order, which is derived from 
equation (ED. From the value of H - A — fn ) / Aq •>the value of cx^m 

satisfying equation ED is obtained. 


we can formally reconstruct the radial velocity as a continuous 
function of time from photometry alone: 


ttrad,9Jt it) 


flc 

Do 


SOI 

ma^m cos(mfif + d m ) 

m=1 


Q.c 

ujo 


m 

yy ma£m cos 

171=1 


mn (t + ^) 


+ (dm 


md\) 


,(45) 


where DJI is the number of components in the multiplet. In the lim¬ 
iting case of DJI —^ oo, this is the time-varying radial velocity com¬ 
ponent, added to the systemic velocity, that would be measured by 
spectroscopy. The number of terms, DJI, is, however, very limited 
and not large enough in practical cases, hence we need to devise a 
more practical way of solving the problem. Nevertheless, equation 
is instructive in the sense that it proves that the radial velocity 
can be reproduced from the frequency spectrum of the photometric 
variations. 


4 FM METHOD 

The amplitudes and phases of the multiplet components are deter¬ 
mined by the binary orbit. Hence, inversely, once we have observed 
the frequency spectrum of a pulsating star in a binary system, we 
can extract infor mation of the binary orbi t. Essentially, this is the 
the FM method dshibahashi & Kurtal2012l) . 

While the strict values of the orbital elements can be obtained 
by analysing the continuous radial velocity, [equation 145 H. if the 
number of components of the multiplet DJI is large enough, in prac¬ 
tice DJI is limited to 3 or 4. With such a limited number of terms, 
equation G5} may still substantially deviate from the real radial ve¬ 
locity. In order to improve this, we derive the orbital parameters 
first and reconstruct the radial velocity in the case of DJI —» oo. 
It is then instructive to derive some ‘easy solutions' of the orbital 
parameters to illustrate the technique. More approximate solutions 
can be obtained directly from the Fourier transform of the light 
curve without carrying out further numerical computations. In this 
section, we treat such cases. 


Figure 5. TZrn (a , e, zu). for a fixed value of a = 0.217, as a function of 
e. The width of each curve shows the range of zu from 0 to 271. 


4.1 The eccentricity, e 

Equation 442b leads to the following relation between the ampli¬ 
tudes of the m + 1-th and m-th sidelobes: 

A+(m+l) + + _ Jl« m +l) 

A — m ~\~ A— m Jo (rr^m+1) 

=\Hm(a,e,vj). (47) 

The RHS of equation (1461 . defined as 7Z(a, e, zu), is a function of 
a, e and zu, while the LHS is observable. Fig.[5]shows 7 Z(a, e, zu) 
as a function of e, for m = 1, 2, and 3, for a fixed value of a = 
0.217. The width of each curve shows the range of zu from 0 to 27t. 


' JiK„ 

Jo(ot£.r, 


(46) 


In the case of af m < 1, Ji(a£m) — a^m/2 and Jo(a£ m ) ~ 
1, hence 

Um ~ (48) 

S m 

and then TZ is regarded as being no longer dependent on a. It should 
be noted that, while {£m} are dependent on e and zu, the ratios 
{(.m+i/fm}, and then TZ m , are dependent mainly on e, but rela¬ 
tively weakly on zu (see Fig.[5j. A good initial guess for the eccen¬ 
tricity is then derived from the amplitude ratio of the (m + l)-th 
sidelobes and m-th sidelobes. Indeed, as seen in equation ( fl4t . the 
ratio £ m +i/£m is, in a good approximation, given by 

£m+1 _ m Jm+i((m + l)e) 

fm — m + 1 Jm(me) 

Whether or not a£ m <g 1 can be checked by seeing the amplitude 
ratio of the sidelobe to the central component: 


A+m A — m _ 

Ao Jo(a£m) 

= : Sm(o/.^m). 


(50) 

(51) 


If < 1, Sm(at)m) — otfm <C 1. So, if the ampli¬ 

tude ratio is substantially smaller than unity, the assumption of 
ctf m < 1 is justified. In such a case, the eccentricity is estimated 
as dShibahashi & Kurt3l2012l) 


A+2 + A -2 

6 7^ 2 — j • 

A+i+A-z 


(52) 


Use of equation 449b for more than two values of m deter¬ 
mines multiple values of e, each with its own uncertainty. The final 
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value of e should be determined as a weighted average of those, 
which, given that A m decreases as m increases, will lead to the 
eccentricity being dominated by the £1 and £2 terms. 

4.2 The angle between the periapsis and the nodal point, vj 

The angles $1 and $2 are regarded as functions of w alone, since 
the eccentricity e has already been obtained. The angle zu between 
the periapsis and the nodal point is then derived from equation GD- 
In the case of e 1, as seen in Fig. 2, the e-dependence of 
both i?i and $2 is weak, and 2$i — $2 — zu. Hence, in such a case, 

w - [0+2(*o) - 4>-2(to)] + (53) 

The uncertainty in 0+2(to) — 0-2(to) by 2n7t can be eliminated 
by taking the inequality of A±i and the deviation of 0±i(to) — 0o 
from 7t/2 into account. 

It should be noted that the angle zu is obtained from infor¬ 
mation of the second sidelobes. If the amplitudes of A ±2 are not 
statistically significant, we have to conclude that the orbit is close 
to a circular one. 



frequency d 


Figure 6. An amplitude spectrum for the Q5 SC data of KIC 9651065. The 
four highest amplitude peaks were studied in LC. There are no significant 
peaks above 50 d —1 . 

e^ 1 -* by adopting £( 1 \ and then zu^\ We repeat these processes un¬ 
til we get consistent solutions. Fig.[3]is not used in this process, but 
it is illustrative. 


4.3 The projected semi-major axis, a± sin i 

By substituting the values of e and zu thus obtained into the first 
line of equation CED. we estimate £1 (e, zu). The value of a is then 
obtained by dividing (equation 15 1H by £ m . Then, the pro¬ 
jected semi-major axis is derived by equation 420l >: 

CL 

a\ sini = —c. (54) 

+!q 


4.4 The mass function, /(mi, m 2 , sin i) 


Once the value of ai sin i is derived, with the help of the pulsation 
frequency and the orbital period, the mass function can be derived: 


/(mi, m 2 , sini) := 


m -2 sin 3 i 
(mi + m2) 2 


^rb 27TG’ 


(55) 


where G is the gravitational constant. With a suitable asteroseis- 
mic estimate (or a reasonable assumption) of the primary mass, the 
minimum mass of the secondary star, m 2 sin i, can be deduced. 


5 EXAMPLES WITH REAL KEPLER DATA 

The Kepler mission accumulated time series of photometric data 
of about 200 000 stars from 2009 — 2013. The long time base and 
high duty cycle of Kepler data make them the data of choice for 
FM analyses. Here we present examples for three highly eccentric 
Kepler targets: KIC 9651065, KIC 10990452 and KIC 8264492. In 
each case, we use multi-scale, maximum a posteriori (msMAP) 
pipeline data collected with long cadence (LC, 29.45-min) sam¬ 
pling. 


5.1 KIC 9651065 

KIC 9651065 is a K v = 11.1 binary system showing first, second 
and third FM sidelobes, hence it has high eccentricity. According to 
iHuber et~ai](l2014l) . its effective temperature and gravity are T e g = 
7014+ 1 g| K and logg = 3.8331 q JJ 2 (cgs units), and it mass is 
mi = 1 . 703 ±ons 5 M 0 , which we round to T efi = 7400 ± 150 K, 
log g = 3.83 ±0.13 and mi = 1.70±0.17M o 


4.5 The radial velocity curve deduced from photometry 

Once e and zu are determined, we can deduce cos / and sin / 
from equations CQ} and G3- respectively, as functions of the mean 
anomaly, ff(f — f p ), all the terms of £n(e, zu) and d(e, zu). Then 
substituting them, along with a sin i and H, into equation |7J, we 
obtain the radial velocity as a function of the orbital phase. 

4.6 Iteration for an improvement of parameters 

Once the value of zu is estimated, we may iterate the above men¬ 
tioned processes, by taking account of the tii-dependences of £ m . 
The value of zu newly derived in this way may be slightly differ¬ 
ent from that obtained previously. The iteration should be repeated 
until consistent solutions are obtained. An initial guess for £ is ob¬ 
tained by ignoring ^-dependence: = £(e^°\ zu = 0). The first 

iteration gives = ^(e^°\zu^). We improve the value of e as 


5.1.1 Kepler mission data 

The Kepler mission data available for this star are QO-17 LC data, 
and Q5.1-5.3, Q7.1-7.3 short cadence (SC) data with 1-min inte¬ 
gration time. An examination of the Q5 SC data shows that there 
are no significant peaks with amplitudes greater than 10 pniag at 
frequencies higher than 60 d" 1 . In the frequency range 0 — 50 d" 1 
the four highest amplitude peaks were studied for FM. Fig.[6]shows 
the amplitude spectrum for the Q5 SC data in the frequency range 
0 — 50d _1 where the highest peaks chosen for analysis can be 
seen. They are in order of amplitude: v\ = 19.47768 d -1 , V 2 = 
30.80189 d" 1 , u 3 = 21.71214 d' 1 and u A = 36.14643 d" 1 . 
There is no clear separation of pmode and gmode frequency 
ranges; there are peaks in all frequency ranges, suggesting p modes, 
g modes and mixed modes. 

The Q0, Q1 and Q17 LC data were of shorter duration than full 
quarters, and they had a significant slope from the pipeline reduc¬ 
tions. We therefore chose to analyse the Q2-16 LC data covering 
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Table 1. A least-squares fit of the frequency septuplet for the highest amplitude mode to the Q2-Q16 Kepler data for KIC 9651065. The frequencies of the 
multiplet are separated by the orbital frequency, v orb = 0.0036524 ± 0.0000036 d — 1 (P OT h = 273.8 it 0.3 d). Column 4 shows that the first sidelobe phases 
are close to rr/2 = 1.57 rad out of phase with the central peak, as required by theory for FM. The zero point for the phases has been chosen to be a time when 
the phases of the first sidelobes are equal, to = B JD 2455783.05262, for the first multiplet. Column 5 shows that the phases of the first sidelobes of the other 
multiplets are equal within the errors at this time. 


frequency 

d” 1 

amplitude 

mmag 

phase 

rad 

(0± 1 ) - 00 

rad 

4 > +m P—m 

rad 

A| m +A m 

A +rn -A rn 

A o 

A_|_ m +A_ m 

19.46672 

19.47037 

19.47402 

19.47768 

19.48133 

19.48498 

19.48863 

0.0223 ± 0.0020 
0.0613 ± 0.0020 
0.2213 ±0.0020 
1.9308 ± 0.0020 
0.2153 ±0.0020 
0.0443 ± 0.0020 
0.0152 ±0.0020 

-1.4619 ±0.0896 
0.7396 ±0.0326 
-3.1416 ±0.0090 
-1.5919 ±0.0010 
-3.1415 ±0.0093 
-0.4639 ±0.0452 
2.1830 ±0.1314 

-1.550 ±0.007 

0.000 ±0.013 
-1.204 ±0.056 
3.645 ± 0.132 

0.2261 ± 0.0029 
0.0547 ±0.0015 
0.0194 ±0.0015 

-0.014 ±0.006 
-0.161 ± 0.027 
-0.189 ± 0.077 

21.70118 

21.70483 

21.70848 

21.71214 

21.71579 

21.71944 

21.72309 

0.0096 ± 0.0020 
0.0332 ± 0.0020 
0.1163 ±0.0020 
0.8504 ± 0.0020 
0.1182 ±0.0020 
0.0195 ± 0.0020 
0.0064 ± 0.0020 

0.5094 ± 0.2090 
2.5352 ± 0.0602 
-1.4206 ±0.0172 
0.1291 ±0.0024 

-1.4362 ±0.0169 
1.3598 ±0.1024 
-2.1510 ±0.3116 

-1.558 ± 0.012 

-0.016 ± 0.024 
-1.175 ± 0.119 

-2.660 ±0.312 

0.2758 ± 0.0034 
0.0466 ± 0.0033 
0.0188 ±0.0033 

0.008 ± 0.012 
-0.346 ± 0.076 
-0.200 ±0.180 

30.79094 

30.79459 

30.79824 

30.80189 

30.80555 

30.80920 

30.81285 

0.0184 ±0.0020 
0.0398 ± 0.0020 
0.1291 ±0.0020 
0.7021 ± 0.0020 
0.1365 ±0.0020 
0.0301 ± 0.0020 
0.0075 ± 0.0020 

-1.8457 ±0.1088 
-0.1224 ±0.0503 
2.4607 ±0.0155 
-2.3441 ± 0.0028 
2.3847 ±0.0147 
-0.8970 ± 0.0664 
0.4618 ±0.2682 

-1.516 ± 0.011 

-0.076 ± 0.021 
-0.775 ± 0.083 
2.308 ± 0.269 

0.3783 ± 0.0042 
0.0996 ± 0.0040 
0.0369 ± 0.0040 

0.028 ± 0.011 
-0.139 ± 0.041 
-0.421 ± 0.118 

36.13547 

36.13913 

36.14278 

36.14643 

36.15008 

36.15374 

36.15739 

0.0189 ±0.0020 
0.0389 ± 0.0020 
0.0683 ± 0.0020 
0.2864 ± 0.0020 
0.0659 ± 0.0020 
0.0107 ±0.0020 
0.0084 ± 0.0020 

-2.1439 ±0.1061 
-0.4995 ±0.0515 
2.8955 ± 0.0293 
-1.9325 ±0.0070 
2.8373 ± 0.0304 
-0.7266 ±0.1873 
1.8010 ±0.2372 

-1.484 ±0.022 

-0.058 ± 0.042 
-0.227 ±0.194 
3.945 ± 0.239 

0.4686 ±0.0104 
0.1732 ±0.0099 
0.0953 ± 0.0099 

-0.018 ± 0.021 
-0.569 ± 0.066 
-0.385 ± 0.111 


a time span of 1388.2 d. We begin the discussion of the analysis 
with the highest amplitude peak at 19.47768 d -1 . Fig.[7]shows the 
prewhitening process. The first, second and third FM sidelobes are 
significant and detected. There is some amplitude and/or frequency 
modulation of the v\ mode, hence there is amplitude left near to 
that frequency after prewhitening. This is astrophysical, since am¬ 
plitude and frequency variations on t he time scale of this d ata set 
are common in 8 Set stars (see, e.g jBowman & Kurtzil 20141 ). 

This process was carried out for the other three of the high¬ 
est four peaks. Figures are not shown for these other peaks. The 
results are given in Table[T] In the small sections of the amplitude 
spectrum between significant peaks the highest noise peaks have 
amplitudes of 8 pmag. The error on the amplitude determination 
is taken to be one quarter of that, with phase an d frequenc y er¬ 
rors scaled accordingly iMontgomerv & O’ Donoghueil 19991) . The 
formal least-squares errors contain variance from all the peaks, so 
must be scaled down by this estimate. 

5.1.2 FM analysis 

For all four modes presented in TableQ]all sidelobes are equally 
split from the central frequency, being consistent with theoretical 
expectation of FM stars. This splitting gives the orbital frequency, 
u OTb = 0.0036524 ± 0.0000036 d" 1 (P orb = 273.8 ± 0.3 d), 
where we have taken the best value here to be from the splitting of 


the first sidelobes of the highest amplitude mode. The zero point for 
the phases has been chosen to be a time when the phases of the first 
sidelobes of the highest amplitude peak are equal. The phases of 
the first sidelobes are close to n/2 rad out of phase with the central 
peak. That they are slightly less than this is a consequence of the 
values of the eccentricity and argument of periapsis. The same is 
true for the other modes within the error bars. Furthermore, the 
phase relation (4>+ n — 4 >~ n ) for n = 2, 3, 4 is the same within the 
errors for the u\ and V 2 modes. As for the 1/3 and 1/4 modes, the 
uncertainty is too large to confirm this. These facts are consistent 
with the theory for FM discussed in Section 3.3. 

For each set of sidelobes, the values of ma£ m and 9 m — m0\ 
that were derived by equations ( 142b and d44b . respectively, are sum¬ 
marized in Table[2] The eccentricity e is then estimated from the 
first and second sidelobes around v\ by equation 152} to be e = 
0.569 ± 0.030. The technique explained in Section 4.1 and inspec¬ 
tion of Fig.[3]gives the coefficient £i(e, vj ) = 0.840±0.015. Using 
this value and the amplitude ratio for u\, a£i = 0.2247 ± 0.0015, 
we obtain a = 0.2675 ± 0.0018. The mass function is then 
estimated to be 0.0992 ± 0.0057, which gives a minimum sec¬ 
ondary mass of m 2 = 0.87 ± 0.02 Mq based on an assumption 
of mi = 1.7Mq, so the secondary is probably a main sequence 
G star. We derive the semi-major axis of the primary star about the 
barycentre from equation {20]l to be m sini = 0.378 ± 0.007 au 
from the data for v\. 
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Figure 7. Amplitude spectra for KIC 9651065 for the Q2-16 LC data. The top, left panel shows the highest amplitude pulsation peak at v\ = 19.47768 d —1 . 
The top, right panel shows the amplitude spectrum after prewhitening where two equally spaced sidelobes can be seen easily. These are the first FM 
sidelobes. A careful examination shows the second and third FM sidelobes also. The bottom left panel shows the amplitude spectrum after prewhitening 
and its first FM sidelobes, where the second and third sidelobes can be seen. The remaining variance close to the mode peak after pre-whitening the second and 
third sidelobes is seen in the bottom right panel; this is caused by amplitude modulation of the mode frequency over the time span of the data set, a common 
characteristic of 5 Set stars. The highest remaining peak may be that of an independent pulsation mode. 


Table 2. The estimated quantities of the terms for KIC 9651065, ma£ m 
and 6m, appearing in the expression of the radial velocity. Column 1 shows 
the frequency of the central peak, and Column 2 shows its ratio to the orbital 
frequency. Here c is the speed of light. Column 4 shows raa£ m estimated 
from equation 63- Column 5 shows 6 m — m6\ estimated by equation 63). 


^puls 

d- 1 

(k'orb/^puls)^ 
kms -1 

m 

nia^m 

6 m — m0\ 
rad 

19.47768 

57.0 ± 0.2 

1 

0.2247 ±0.0015 

0.000 ± 0.007 



2 

0.1094 ± 0.0030 

0.969 ± 0.028 



3 

0.0582 ± 0.0045 

4.964 ± 0.066 

21.71214 

51.1 ±0.2 

1 

0.2732 ± 0.0033 

-0.008 ±0.012 



2 

0.0932 ± 0.0066 

0.983 ± 0.059 



3 

0.0564 ± 0.0099 

4.953 ±0.156 

30.80189 

36.1 ±0.1 

1 

0.3717 ±0.0040 

-0.038 ±0.011 



2 

0.1992 ± 0.0080 

1.184 ±0.042 



3 

0.1107 ±0.0120 

4.295 ±0.134 

36.14643 

30.7 ±0.1 

1 

0.4563 ± 0.0096 

-0.029 ±0.021 



2 

0.3452 ± 0.0194 

1.457 ±0.097 



3 

0.2856 ± 0.0297 

5.114 ±0.120 


The facts that A +2 — A -2 < 0 and (</>±i) — (fro > — n/2 
indicate 0 < 2i?i — #2 < 7t. The phase difference between 
(j > + 2 and (j>- 2 , along with this constraint, leads to 2$i — $2 = 
2.17 ±0.03 rad. With use of the derived value of e, this gives 
t*7 = 2.22 ± 0.04 rad. 

Fig.[8]shows the radial velocity curve of KIC 9651065 derived 
from the Kepler photometric data alone. This demonstrates again 


that the present FM method extracts binary information from the 
Kepler light curve without spectroscopic observations. 


5.1.3 Comparison with the result obtained by PM 

iMurphv et alj d2014l ) developed a different method, by which the 
phase modulation (PM) in the time domain of intrinsic pulsation 
frequencies of the star is tracked. In this method, first of all, the 
frequency of the central component of the multiplet is measured 
in the Fourier transform. Then the light curve is divided into short 
segments and, with the frequency fixed, the pulsation phase in each 
segment is measured by a least-squares method. This provides us 
with “time delays” r(t) as a function of time, where r(t) is defined 
as in equation 0. Then by carrying out a Fourier analysis, we ob¬ 
tain the amplitudes au!g 1 ^n(e, tzr) and phases d n (e,zv) given in 
equation l !18t for n = 1, 2, 3, • ■ The amplitude ratios between £ n 
and £n+i allow us to estimate the eccentricity e. On the other hand, 
the radial velocity is deduced as a function of time by taking the 
time derivative of the time delay. 


Vrad (t) = c—. (56) 

Then, from the maximum and minimum radial velocity, along with 
the eccentricity deduced from the amplitude ratios, the projected 
semi-major axis is deduced; 

til Sin i — \/1 C^(^rad,max tt ra d,min)' (57) 

The angle 7*7 is also deduced from w ra d,max and u ra d,min: 
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Orbital phase 


Figure 8. The radial velocity curve, derived from the Kepler photometric 
data with the FM method, of KIC 9651065. The zero point for the orbital 
phase corresponds to the periapsis passage of the star. 



frequency d 


Figure 9. Amplitude spectrum of the LC Ql-16 data for KIC 10990452. 
Most of the variance is in the pmode frequency range with a dominant 
highest peak at v\ = 17.72374 d —1 . 


Table 3. Comparison of the binary parameters of KIC 9651065 derived with 
the FM method and the PM method. 


Quantity 

PM 

FM 

Porb M) 

272.70 ±0.82 

273.80 ±0.30 

e 

0.47 ±0.03 

0.57 ±0.03 

zu (rad) 

2.01 ±0.30 

2.22 ±0.04 

ai sin i (au) 

0.37 ± 0.02 

0.38±0.01 

/(mi, m 2 , sin i) (Mq) 

0.0916 ±0.0108 

0.0992 ±0.0057 


1 ^rad.max H - ad,min /en\ 

COS VJ = --- (58) 

6 L^radjmax ^rad,min_ 

iMurphv et alj i2014l ) applied this PM method to sev¬ 
eral pulsating stars observed by the Kepler Mission, including 
KIC 9651065. They derived the orbital elements of this star, as well 
as the radial velocity curve. The results obtained in this paper by 
the FM method and their PM results are in good agreement (see 
Tabled). 

It shoul d be noted here that the PM m ethod has been further 
developed bv lMurphv & Shibahashil d2015ti . The latest method de¬ 
rives directly the orbital elements without converting time delays to 
radial velocities, offering higher precision. 

5.2 KIC 10990452 

KIC 10990452 is a K p = 12.4 binary system showing first, second, 
third and one component of the fourth FM sidelobes, indicatin g tha t 
the system is highly eccentric. According to iHuber et al.l d2014h . 
its effective temperature and gravity are T e g = 75851287 K 
and log 5 = 4.0221°;^® (cgs units), and its mass is m = 
1 . 691 !q 219 Mq, which we round to T e g = 7600±250 K, logs = 
4.02 ± 0.13 and mi = 1.69 ± 0.33 Mq. These are consistent with 
the original KIC parameters. We used Ql-16 data in our analy¬ 
sis, which have gaps because the star fell on the failed Module 3. 
Light curves in quarters Q7, Qll and Q15 are therefore missing. 
The short 9.7-d engineering Q0 and the short 31-d final Q17 were 
not included. 

Fig-ED shows the LC amplitude spectrum where it can be seen 
that there is a dominant highest peak. We examine that and the sec¬ 
ond highest amplitude peak next to it for FM. The top left panel 
of Fig.flOl shows the highest peak at v\ = 17.72374 d' 1 and the 
second highest peak close to it, vi = 17.85683 d' 1 . Because of 


the complexity of the spectral window caused by the gaps in the 
data, the window patterns for the two frequencies interfere, so they 
must be modelled together. The top right panel shows the amplitude 
spectrum after the highest two mode peaks have been prewhitened; 
the FM sidelobes are apparent for both modes. The middle left 
panel is one step further with prewhitening of the first sidelobes, 
too. The remaining central amplitude of V 2 represents amplitude 
modulation over the time span of the data set - either astrophys- 
ical or instrumental. Given that v\ does not show this as much, 
it is probably astrophysical. The middle right panel shows the third 
sidelobes and the low-frequency component of the fourth sidelobes. 
The bottom panel shows the amplitude spectrum of the residuals af¬ 
ter the two nonuplets given in Table[4]have been prewhitened. 

When searching for significant peaks associated with pulsa¬ 
tion modes among thousands of independent frequencies in the am¬ 
plitude spectrum - as we do in this case - it is found that the highest 
noise peaks have amplitudes about 4 times that of the rms noise in 
amplitude. The least-squares fitting that we do to find frequencies, 
amplitudes and phases estimates the errors based on the total vari¬ 
ance in the data. Since there are so many other pulsation frequen¬ 
cies present that have not been modelled, the noise is not white and 
errors are significantly overestimated. We therefore looked at a sec¬ 
tion of the amplitude spectrum that is pure noise and found that the 
highest noise peaks are around 12 pmag. The amplitude error was 
therefore scaled to one quarter of that, or 3 prnag. The phase and 
frequency errors were scaled by the same factor, since the errors on 
those quant ities are proportion al to the amplitude signal-to-noise 
ratio ^Montgomery & O' Donoghudl 19991) . 

The results of the least-squares fits of the two pulsation mode 
peaks and their first, second, third, and fourth FM sidelobes are 
given in Table[4] All sidelobes are equally split from the central 
frequency, and is consistent with theoretical expectation for FM 
stars. This splitting gives the orbital frequency, v ar b = 0.008190± 
0.000024 d _1 (Porb = 122.11 ± 0.36 d). The zero point for the 
phases has been chosen to be a time when the phases of the first 
sidelobes of the highest amplitude peak are equal. The phases of 
the first sidelobes are 7t/2rad out of phase with the central peak. 
The same is true for the second highest amplitude mode with V 2 - 
Furthermore, the phase relation (0+ m — <t>-m ) for m = 2, 3, 4 
is the same within the errors for the two modes. These facts are 
consistent with theory for FM discussed in Section 3.3. 

For each set of sidelobes, the values of and 9 m — m9i 

were derived by equations §2) and 444t . respectively, and are sum- 
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Figure 10. Top left: Amplitude spectrum for KIC 10990452 in the frequency range of the two highest amplitude p-mode peaks. The spectral window has 
significant sidelobes because of the missing quarters of data. Top right: This panel shows the first FM sidelobes after prewhitening by and V 2 ; further 
sidelobes can be seen. Bottom left: The this panel shows the data after prewhitening by the first sidelobes, too, where the second and third sidelobes are easily 
visible. Bottom right: In this panel the third sidelobes are easily seen for both main frequencies, and the low frequency fourth sidelobe is also seen for the ui 
group. 


Table 4. A least-squares fit of the frequency nonuplets for the two highest amplitude modes to the Q1 to Q16 Kepler data for KIC 10990452. The frequencies 
of the multiplet are separated by the orbital frequency, is 0 rb = 0.008190 d= 0.000024 d -1 (P Q rb = 122.11 d= 0.36 d). The zero point for the phases has been 
chosen to be a time when the phases of the first sidelobes for the highest amplitude peak are equal, to = B JD 2455673.01863. Column 4 shows that the first 
sidelobe phases are n/2 = 1.57 rad out of phase with the central peak, as required by theory for FM. Column 5 shows that the phases of the first sidelobes are 
equal within the errors at this time. 


frequency 

amplitude 

phase 

(<t>± l) ~ 00 

m 

Aj rrn -\-A_ rn 

A +m -A_ rn 

Aq 

A-\-rn~^~ A — m 

d" 1 

mmag 

rad 

rad 

rad 


17.6909790 

0.018 ± 0.003 

0.7772 ±0.1691 





17.6991685 

0.021 ± 0.003 

-0.1864 ±0.1450 





17.7073581 

0.044 ± 0.003 

-1.0048 ±0.0684 





17.7155476 

0.194 ±0.003 

-2.1688 ±0.0154 





17.7237371 

5.654 ± 0.003 

-0.6296 ± 0.0005 

-1.539 ±0.012 




17.7319266 

0.168 ± 0.003 

-2.1688 ±0.0179 


0.000 ± 0.024 

0.064 ± 0.001 

-0.072 ±0.012 

17.7401161 

0.052 ± 0.003 

3.1176 ±0.0577 


4.122 ±0.090 

0.017 ±0.001 

0.083 ± 0.044 

17.7483056 

0.027 ± 0.003 

1.6975 ±0.1114 


1.884 ±0.113 

0.008 ± 0.001 

0.125 ±0.089 

17.7564952 

0.010 ±0.003 

0.4999 ±0.3118 


-0.277 ±0.312 

0.005 ± 0.001 

-0.286 ±0.158 

17.8240752 

0.002 ± 0.003 

-2.5761 ± 1.3006 





17.8322647 

0.011 ± 0.003 

2.1708 ±0.2699 





17.8404542 

0.026 ± 0.003 

1.6322 ±0.1155 





17.8486437 

0.095 ± 0.003 

0.4765 ± 0.0317 





17.8568332 

2.674 ± 0.003 

1.9836 ±0.0011 

-1.501 ±0.024 




17.8650228 

0.083 ± 0.003 

0.4883 ± 0.0362 


0.012 ±0.048 

0.067 ±0.002 

-0.067 ±0.024 

17.8732123 

0.029 ± 0.003 

-0.6460 ±0.1022 


-2.278 ±0.154 

0.021 ±0.002 

0.055 ± 0.077 

17.8814018 

0.009 ± 0.003 

-1.8700 ±0.3233 


-4.041 ± 0.325 

0.007 ±0.002 

-0.100 ±0.213 

17.8895913 

0.004 ± 0.003 

-2.4669 ± 0.6935 


0.109 ±0.694 

0.002 ± 0.002 

0.333 ± 0.745 
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Orbital phase 

Figure 11. The radial velocity curve, derived from the Kepler photometric 
data with the FM method, of KIC 10990452. The zero point for the orbital 
phase corresponds to the periapsis passage of the star. 

marized in Table[5] The values of these quantities derived from u\ 
and V 2 are in good agreement, and this is consistent with theoret¬ 
ical expectation of FM stars. The eccentricity e is then estimated 
from the first and second sidelobes around by equation {Uto 
be e = 0.569 ± 0.030. 

The technique explained in Section 4.1 and inspection of 
Fig.QJgives the coefficient £i(e,-n7) = 0.815 ± 0.015. Using this 
value and the amplitude ratio for v\, cr£i = 0.0640 ± 0.0015, 
we obtain a = 0.0785 ± 0.002. The mass function is determined 
from these values to be / (mi, m 2 , sini ) = 0.0163 ± 0.0011 Mq. 
Assuming mi = 1.7 Mq llHuber etah l2014h gives a minimum sec¬ 
ondary mass of m 2 = 0.41 ± 0.01 Mq, so the secondary is prob¬ 
ably a main sequence M or K star. We also derive the semi-major 
axis of the primary star about the barycentre. From equation 120b . 
we find ai sin i = 0.122 ± 0.003 au from the data for v\. 

The facts that A+ 1 — A- 1 < 0 and (</±i) — <j >0 < —7t/2 in¬ 
dicate that 37t/2 < 2t?i — t ?2 < 27t. The phase difference between 
rj )+2 and rj>- 2 , along with this constraint, leads to 2i?i — $2 = 
5.79 ± 0.05 rad. With use of the derived value of e, this gives 
w = 5.85 ± 0.05 rad. 

Fig-Qj] shows the radial velocity curve of KIC 10990452 de¬ 
rived from the Kepler photometric data alone. This demonstrates 
that the present FM method extracts binary information from the 
Kepler light curve without spectroscopic observations. 

5.3 The highly eccentric star KIC 8264492 

5. 3.1 PM analysis 

KIC 8264492 is a binary with a <5 Set star component. iHuber et~ai] 
i2014l) give T eff - 7992±3?ji K, logg = 3.947±oi24 (cgs units), 
and M = 1.87 Mq . The bi nary nature of KIC 8 264492 was discov¬ 
ered using the PM method (Mur phy et al.l2014h . The Fourier trans¬ 
form of its light curve (Fig.l 1 2h ) shows a dense frequency spec¬ 
trum, in which the highest amplitude peak occurs at 31.29 d -1 . 
This peak is distinguished from its Nyquist al ias a t 17.7d/ 1 _by 
the larger amplitude of the former iMurphv, Shibahashi & Kurta 
12013b . The non-sinusoidal nature of the time delays indicated the 
orbit is highly eccentric, which is co nfirmed by the high amplitud e 
of the harmonics (see figure 21 of [Murphy & Shibahashi! 120150 . 
From the amplitude ratios of those harmonics, a graphic solution 
for the eccentricity, e ~ 0.73 ± 0.05, is obtained with the help of 
Fig-GD 


Table 5. The estimated quantities of the terms for KIC 10990452, 
and dm, appearing in the expression of the radial velocity. They are esti¬ 
mated from the frequency nonuplets for the two highest amplitude modes, 
u\ and v 2 , shown in Table[4] Column 1 shows the frequency of the cen¬ 
tral peak, and Column 2 shows its ratio to the orbital frequency. Here c is 
the speed of light. Column 4 shows ma£ m estimated from equation <42}- 
Column 5 shows 6m estimated by equation El- 


^puls 

d- 1 

( ^orb / ^puls ) C 
kms -1 

m 

ma£m 

dm — m6l 
rad 

17.7237371 

138.5 ±0.4 

1 

0.064 ±0.001 

0.000 ± 0.012 



2 

0.034 ±0.002 

3.632 ± 0.045 



3 

0.025 ± 0.003 

4.084 ± 0.056 



4 

0.020 ± 0.004 

4.574 ±0.156 

17.8568332 

137.4 ±0.4 

1 

0.067 ±0.002 

0.006 ± 0.024 



2 

0.042 ± 0.004 

3.573 ±0.077 



3 

0.021 ±0.006 

4.263 ± 0.162 



4 

0.008 ± 0.008 

4.767 ± 0.347 


Table 6. The binary parameters of KIC 10990452 derived with the FM 
method. 


Quantity 

PM 

FM 

P OI b (d) 

122.10 ± 0.21 

122.11 ±0.36 

e 

0.55 ±0.03 

0.57 ± 0.03 

vd (rad) 

5.81 ±0.05 

5.85 ±0.05 

a\ sin i (au) 

0.123 ±0.016 

0.122 ±0.023 


5.3.2 FM analysis 

The following FM analysis is based on the highest amplitude peak 
front Fig. EUi. We removed just over 5000 points from the start of 
the Q0-Q16 LC data set because of some trends and flux discon¬ 
tinuities in the early quarters. The data set we used has 64031 LC 
points, and is 1318.9-d long (BJD 2455 072.035 to 2456 390.959). 

Fig-IUJi zooms in on the highest peak at 
v i = 31.2920323 d -1 . The first sidelobes are clearly visible, 
separated from the central peak by 0.004 d -1 . Fig.ll3b shows the 
amplitude spectrum after the central peak (at the dashed red line) 
has been prewhitened. The FM sidelobes are apparent, as shown 
by arrows, at exact multiples of the orbital frequency distant from 
the central peak. 

Table [7] gives the frequencies, amplitudes, phases and their 
uncertainties for the three highest amplitude modes and their de¬ 
tectable pairs of sidelobes. We determined the orbital frequency 
splitting from the splitting of the first pair of sidelobes with respect 
to the central component, u OI b = 0.0039621 ± 0.0000087 d -1 , 
giving P or b = 252.39 ± 0.56 d. 

For each set of sidelobes, the values of ma£ m and 9 m — m6\ 
were derived by equations and m, respectively, and are 
summarized in Table[8] The eccentricity e is then estimated from 
the first and second sidelobes around v\ by equation <123 to be 
e = 0.761 ± 0.045. 

The technique explained in Section 4.1 and inspection of 
Fig.[3] gives the coefficient £i(e, zu) = 0.630 ± 0.030. Using 
this value and the amplitude ratio for vi, afi = 0.369 ± 0.004, 
we obtain a = 0.577 ± 0.028. The mass function is determined 
from these values to be /(mi, m 2 , sin i) = 0.274 ± 0.040 Mq. 
Assuming mi = 1.87 Mq gives a minimum secondary mass of 
m 2 = 1.44 ± 0.10 Mq. We also derive the semi-major axis of 
the primary star about the barycentre. From equation d20b . we find 
a\ sin i = 0.508 ± 0.024 au from the data for v\. 
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Frequency (d' 1 ) 

Figure 12. Fourier transform for KIC 8264492 from 5.0 to 43.9 d —1 . The highest peak is real and is used in our FM analysis. The frequencies used for the PM 
analysis are v\ ... ug = 31.29, 23.92, 33.76, 34.21, 23.54, 21.14, 24.47, 20.99 and 13.53 d 1 . The dashed red line represents the Nyquist frequency of the 
Kepler LC data. 


Table 7. A least-squares fit of the frequency undecuplet for the three highest amplitude modes in the KIC 8264492 data set. The zero-point in time was set 
at BJD = 2455 754.1266, so that the first sidelobes of the highest amplitude mode had exactly equal phases. From the first pair of sidelobes, we obtain a 
frequency splitting of 0.0039621 zb 0.0000087 d 1 , giving P Q rb = 252.39 zb 0.56 d. Frequency uncertainties were calculated from a separate non-linear 
least-squares calculation where sidelobes were not forced to be equally split by the orbital frequency, and should therefore be taken as representative, only. 
Since only three modes were investigated, there is substantial variance left in the data. Treating the remaining variance would reduce the errors. The fact that 
the differences between the phases of the first sidelobes and the central peak for each mode are so close to rt/2 confirms the binarity within the FM theory. 


m 

Frequency 

Amplitude 

Phase 

(<t>± l) - '/’O 

4 > -\-m $—m 

Aj rrn -\-A rn 

A +rn -A rn 

A o 

A +m+ A —m 


a- 1 

mmag 

rad 

rad 

rad 

-5 

31.2722205 ± 0.0001598 

0.009 ± 0.003 

1.495 ±0.380 





-4 

31.2761829 ± 0.0001108 

0.013 ±0.003 

1.065 ±0.263 





-3 

31.2801453 ± 0.0000539 

0.027 ± 0.003 

0.244 ±0.128 





-2 

31.2841076 ± 0.0000300 

0.048 ± 0.003 

-0.115 ±0.071 





-1 

31.2880700 ± 0.0000086 

0.167 ±0.003 

-0.239 ±0.020 





0 

31.2920323 ± 0.0000016 

0.893 ± 0.003 

1.336 ±0.004 

-1.575 ±0.015 




+1 

31.2959947 ± 0.0000085 

0.168 ±0.003 

-0.239 ±0.020 


0.000 ± 0.028 

0.375 ± 0.005 

0.003 ±0.013 

+2 

31.2999571 ± 0.0000229 

0.063 ± 0.003 

-1.083 ±0.054 


-0.968 ±0.089 

0.124 ±0.005 

0.135 ±0.039 

+3 

31.3039194 ± 0.0000330 

0.043 ± 0.003 

-2.138 ±0.078 


-2.382 ±0.081 

0.078 ± 0.005 

0.229 ±0.062 

+4 

31.3078818 ± 0.0000984 

0.015 ±0.003 

-2.579 ± 0.232 


-3.644 ± 0.232 

0.031 ± 0.005 

0.071 ±0.152 

+5 

31.3118441 ± 0.0001587 

0.009 ± 0.003 

-2.401 ± 0.375 


-3.896 ± 0.376 

0.020 ± 0.005 

0.000 ± 0.236 

-2 

23.9104185 ± 0.0000399 

0.036 ± 0.003 

-2.877 ±0.094 





-1 

23.9143809 ± 0.0000144 

0.099 ± 0.003 

2.654 ± 0.034 





0 

23.9183432 ± 0.0000022 

0.656 ± 0.003 

-2.102 ±0.005 

-1.624 ±0.025 




+1 

23.9223056 ± 0.0000150 

0.095 ± 0.003 

2.461 ± 0.036 


-0.193 ±0.050 

0.296 ±0.007 

-0.021 ±0.022 

+2 

23.9262680 ± 0.0000441 

0.032 ± 0.003 

1.861 ±0.105 


4.738 ±0.141 

0.104 ±0.006 

-0.059 ±0.062 

-4 

33.7448167 ±0.0001495 

0.010 ±0.003 

-2.397 ±0.353 





-3 

33.7487790 ±0.0001168 

0.012 ±0.003 

-2.600 ±0.275 





-2 

33.7527414 ± 0.0000533 

0.027 ± 0.003 

-2.810 ±0.125 





-1 

33.7567037 ± 0.0000156 

0.092 ± 0.003 

-2.975 ±0.037 





0 

33.7606661 ± 0.0000034 

0.426 ±0.003 

-1.401 ±0.008 

-1.625 ±0.028 




+1 

33.7646285 ± 0.0000163 

0.088 ± 0.003 

-3.076 ± 0.039 


-0.101 ±0.054 

0.423 ±0.010 

-0.022 ±0.024 

+2 

33.7685908 ± 0.0000468 

0.031 ± 0.003 

2.520 ±0.110 


5.330 ±0.167 

0.136 ±0.010 

0.069 ± 0.073 

+3 

33.7725532 ± 0.0001162 

0.012 ±0.003 

1.799 ±0.277 


4.399 ±0.280 

0.056 ±0.010 

0.000 ±0.177 

+4 

33.7765155 ± 0.0002028 

0.007 ± 0.003 

0.969 ± 0.479 


3.396 ± 0.479 

0.040 ±0.010 

-0.176 ±0.253 
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Figure 13. (a): The amplitude spectrum for KIC 8264492, centred and zoomed on the highest peak. Sidelobes are visible without prewhitening the central 
peak, but in panel (b) the same region is shown with the central peak prewhitened (at the dashed red line). Arrows are drawn at exact multiples of the orbital 
frequency distant from the central peak, showing that four pairs of sidelobes are easily identified. 


Table 8. The estimated quantities of the terms for KIC 8264492, raa£ m 
and Gm, appearing in the expression of the radial velocity. They are esti¬ 
mated from the frequency nonuplets for the three highest amplitude modes, 
v±, V 2 and 1 / 3 , shown in Table[7] Column 1 shows the frequency of the cen¬ 
tral peak, and Column 2 shows its ratio to the orbital frequency. Here c is 
the speed of light. Column 4 shows ma£ m estimated from equation EJ. 
Column 5 shows 6m estimated by equation J44). 


Spills 

d- 1 

(^orb/^puls)^ 

kms -1 

m 

mot <£ 

m 

dm — mOi 
rad 

31.2920323 

38.0 

± 

0.1 

1 

0.369 

± 

0.004 

0.000 

± 

0.014 





2 

0.248 

± 

0.010 

1.087 

± 

0.045 





3 

0.234 

± 

0.015 

1.951 

± 

0.040 





4 

0.124 

± 

0.020 

2.890 

± 

0.116 





5 

0.100 

± 

0.025 

4.335 

± 

0.188 

23.9183432 

49.7 

± 

0.1 

1 

0.293 

± 

0.007 

-0.097 

± 

0.025 





2 

0.208 

± 

0.012 

3.940 

± 

0.071 

33.7606661 

35.2 

± 

0.1 

1 

0.414 

± 

0.009 

-0.051 

± 

0.027 





2 

0.272 

± 

0.020 

4.236 

± 

0.083 





3 

0.168 

± 

0.030 

5.344 

± 

0.014 





4 

0.160 

± 

0.040 

6.395 

± 

0.240 


The minimum mass derived above puts the star in the mid- 
F range, or more massive. That suggests that some of the peaks 
seen in Fig.ll2h may come from pulsations in the secondary. Flow- 
ever, no evidence of anti-phase time delay variations in p modes are 
found, and g modes are low amplitude and not PM sensitive. 

The facts that A +2 —A -2 > Ofor^i and (4>±i) — (j>o < — 7t/2 
for V 2 indicate 37t/2 < 2$i — $2 < 27t. The phase difference 
between (f>+2 and <j>-2, along with this constraint, leads to 2t?i — 
t ?2 = 5.20 ± 0.04 rad. With use of the derived value of e, this gives 
vj = 5.28 ± 0.04 rad. 

The upper panel of Fig.|T4] shows the radial velocity curve 
of KIC 8264492 derived with the FM method from the Kepler 
photometric data alone. The time delay is also derived with the 
FM method. The lower panel of Fig.|T4] shows the variation in 
the light arrival time with respect to the barycentre. This is in 
good agreement with time d elays derived by the PM method (see 
[Murphy & Shibahashil J2015h ). 



Orbital phase 



Orbital phase 

Figure 14. Upper panel: The radial velocity curve, derived from the Ke¬ 
pler photometric data with the FM method, of KIC 8264492. The zero point 
for the orbital phase corresponds to the periapsis passage of the star. Lower 
panel: The time delay of KIC 8264492, derived with the FM method. The 
time delay shown here is the difference in the light arrival time, compared 
with the case where the star lies at the barycentre of the binary system. 

6 DISCUSSION AND CONCLUSIONS 

We have shown in this paper that the Fourier transform of the light 
curve of a pulsating star in a binary system leads to frequency 
multiplets in the amplitude spectra and that all the binary orbital 
information traditionally derived from spectroscopic radial veloc- 
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Table 9. Comparison of the binary parameters of KIC 8264492 derived with 
the FM method and the PM method. 


Quantity 

PM 

FM 

Porb (d) 

253.78 ±1.03 

252.39 ±0.56 

e 

0.67 ±0.04 

0.76 ±0.05 

vj (rad) 

5.28 ±0.05 

5.28 ±0.04 

a± sin i (au) 

0.41 ±0.05 

0.51 ±0.02 

/(mi,m 2 ,sin i) (M 0 ) 

0.143 ±0.054 

0.274 ±0.040 


ity measurements can be derived from the frequency splitting and 
the amplitudes and phases of the components of the frequency 
multiples This is an exte nsion of the FM method described in 
[Shibahashi & Kurt/I 120121 ). in which the case of e ^ 1 was mainly 
discussed. We have improved the applicability to highly eccentric 
cases, and presented a method for determining binary orbital pa¬ 
rameters, including the eccentricity and the argument of periapsis, 
from photometry alone. 

The orbital elements to be determined are (i) the orbital angu¬ 
lar frequency, Q, (ii) the projected semi-major axis, a sin i, (iii) the 
eccentricity, e, and (iv) the argument of the periapsis, vj. Hence, 
we need at least four independent relations among the components 
of the frequency multiplet. The orbital angular frequency is di¬ 
rectly measured from the frequency spacing of the adjacent com¬ 
ponents of the multiplet. The combination of the amplitude ratio 
of the sum of a pair of ± 1-components to the central component, 
(A+i ± A-i)/Aq, and the ratio for the case of ±2-components, 
(A+ 2 ± A- 2 )/Aq, gives ai sin* and e. The phase information of 
the first and the second sidelobes leads to the argument of the pe¬ 
riapsis. Careful examination of the phase differences also indicates 
whether the periapsis is located at the far side or the near side of 
the orbit, with respect to us. Hence, even if both components of the 
binary are pulsating, we can separate the pulsation spectra of two 
pulsating stars. 

Binary information can be extracted by the PM method 
dMurphv et al.l2014l ; lMurphv & Shibahashil2015l ). The present FM 
method and the PM analysis are complementary. The PM analy¬ 
sis offers a clear visualization of the binary orbit in the time do¬ 
main by dividing the light curve into segments, while the present 
FM method utilises the full frequency resolution of the data in the 
Fourier domain, and leads to high precision, particularly in the case 
of short-period binaries. 

For pulsating stars, such as 8 Set stars, with many significant 
peaks in the amplitude spectrum, some care is needed to test for 
possible unresolved contamination of FM multiplet components by 
other, independent mode or combination frequencies. This is easily 
done for multi-periodic pulsators, because the amplitude and phase 
relationships for each FM multiplet must agree within the relations 
given in this paper. In the PM method, unresolved contamination 
results in obvious outliers in the time delay diagram, again making 
it easy recognise and discard contaminated mode frequencies. 

A promising application of our method is in the search for 
exoplanets orbiting pulsating 5 Scuti stars, or more generally, exo¬ 
planets around upper main-sequence stars. Both the transit method 
and the ground-based Doppler method, which are widely used in 
the search for exoplanets around solar-like stars and have discov¬ 
ered most of the known exoplanets to date, are more difficult to 
apply for upper main-sequence stars. Planetary transits are smaller 
in the case of upper main-sequence stars because of the much larger 
stellar disks. On the other hand, the typical rotationally broadened 
spectral lines of A-type and earlier type stars make it harder to de¬ 
tect the Doppler shift of spectral lines. Our new technique of pho¬ 


tometric measurement of binary orbital parameters opens a possi¬ 
bility of exoplanet detection around upper main sequence pulsating 
stars, such as 8 Set stars, and around compact pulsating stars. 

Another application of the present method is at the opposite 
extreme, that is, in the search for invisible massive companions in 
binary systems. If a pulsating 8 Set star forms a binary system with 
a star born initially with mass less than about 8 M@ but more than 
the pulsating star itself, the massive component must have already 
evolved into a white dwarf while the less massive component is 
still in the main-sequence stage as a pulsating 8 Set star. The mass 
of the white dwarf must be less than the Chandrasekhar mass limit, 
but still be of ~ 1 Mg , and the star must be too faint to be detected. 
On the other hand, if a binary system is composed of a <5 Set star 
and a star bom initially more massive than ~ 8 Mg but less than 
~ 30 Mg , the massive star must have already been turned to a neu¬ 
tron star, probably with mass ~ 1 Mg, and the star must be again 
undetectable except perhaps as a pulsar. Although the survival of 
the companion after the explosion of the massive star is uncertain 
from the theoretical viewpoint, we know of binary systems com¬ 
posed of a neutron star and either an early-type massive star, as 
a high-mass X-ray binary (HMXB), or a cool star, as a low-mass 
X-ray binary (LMXB). Hence it seems natural to expect a binary 
system of a neutron star and an A-type 8 Set star. The further ex¬ 
treme case is a star born with an initial mass > 30 Mg. Such a star 
becomes a black hole. A well-known binary system composed of a 
black hole candidate is Cygnus X-l, whose companion is a blue su¬ 
pergiant variable star. If a black hole forms a binary system with a 
8 Set star, the present FM technique is a unique, promising method 
of finding such an exotic system. 
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A.2 Derivation of equation i tTH 

Next, let us consider sin /. From equation © 


sin / = \/l — e 2 — 


ai . 


(A. 11) 


APPENDIX 1: DERIVATION OF EQUATIONS (11) AND (12) 

We need to write the radial velocity as an explicit function of time. 
In doing so, the essential point is how to express trigonometric 
functions of the true anomaly / as an explicit function of time. 
In this Appendix, we derive these expressions. 


A.l Derivation of equation ( HU 

Let us first consider cos /. The distance r between the barycentre 
and the star is expressed with help of a combination of the semi¬ 
major axis ai. the eccentricity e and the true anomaly /: 

ai(l - e 2 ) 


1 + e cos / ’ 

and then 


•' = ;[ (1 - e a )21 - i]. 


(A.l) 


(A.2) 


Since, with the help of the eccentric anomaly u, as derived from 
equations © and ©, 

r = ai(l — ecosw), (A.3) 

the term a\/r in the right-hand-side of equation dA.2b is 


ai _ 1 

r 1 — e cos u 


(A.4) 


Kepler’s equation links the eccentric anomaly u with the time after 
the periapsis passage: 


u — e sin u = f l(t — t p ). 

The partial derivative of u with respect to t leads to 
du ft 

dt 1 — e cos u ’ 


(A.5) 


(A.6) 


and then, with the help of equations dA.4b and dA.6b . equation 1A.21 
is reduced to 


cos / = 


^ e ' ) Qdt 1 


(A.7) 


The eccentric anomaly u is expressed as a solution of Kepler’s 
equation, in terms of Fourier expansion with respect to the time 
after periapsis passage [see t extbooks on celestial mechanics, e.g., 
iBrouwer & Clement ( 1961V : 


o° 1 

, = fi (t — t p ) + 2 — J„(ne) sin nQ(t — t p ). 


Then, 

1 du 
12 dt 


= 1 + 2 Jn(ne) cosnf2(f — t p ). 


(A.8) 


(A.9) 


The partial derivative of u given in equation d A.5I ) with respect to e 
leads to 
du sin u 


de 1 — e cos u 

Then, with the help of equation dA,4b . 
r du 


sin / = a/1 — e 2 —. 

de 

On the other hand, from equation dA.8b . 

du °° 

= 2 ^2 Jn'(ne) sin nQ(t — t p ), 


(A. 12) 


(A. 13) 


(A. 14) 


where J n '(x ) := dJ„(a:)/da;. Hence, we obtain the expression 
given in equation 03: 


sin / = 2-y/1 — e 2 J n '(ne) sinnf2(t — f p ). 


(A. 15) 


APPENDIX 2: A FREQUENCY UNDECUPLET 

In Section [+2l we truncated the infinite series of equation d26b with 
N = 5, and derived the relations between the complex amplitudes 
of the first and second sidelobes. The complex amplitudes of the 
higher order sidelobes of the undecuplet are given as follows. 


+ i 


. f Ji(a^2) \ . 


\Jo(a&) J 


sin(2i9i — $ 2 ) 


(A. 16) 


Aq \Jo(a£ 2 )J 


“Us+yJ v m. 


+ i 


f Ji(a&)\ 1 ( J2(a£i)\ . 


\Jo{a£,2)J \Jo{at;i)J 


sin(2t9r — $ 2 ) 


L17) 


A± 3 _ _j_ / JijaY \ 
Aq \Jo{a>t; 3 )J 


1 ± 


( M^ y 1 

Vo(q6)/ 


( Ji(a£i)Ji(a£2) \ 

V Jo(a£i)Jo(ct£ 2 ) J 


cos ($1 + 1?2 — # 3 ) 


f \ ( 7i(a£i) Ji(q6) \ 
\Jo(u£3) J \Jo(a£i)Jo{a&) J 


x sin(t?i + t >2 — $ 3 ) 


a ± W 3 


(A. 18) 


Substituting equation dA.9b into equation dA,7b . we reach the ex¬ 
pression given in equation £□): 


cos f — —e + —J n (ne) cos nf2(t — f p ). (A. 

p / 


10 ) 


A±4 

Aq 


VJo«4)y 


1 ± 


( JYY \ 1 

V Jo(af 4 )y 


*<(Siy ) c » (2 *-’ ,i) 
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Ji(a£i) Ji(a^3) 
Jo(«£i) Jo{a£ 3 ) 


cos(i?i + 1?3 — $4) 


JiK4)\ ^ 2 ( 0 ^ 2 ) 
Jo ( 0 ^ 4 ) / VJo(a^ 2 ) 


sin(2i?2 — $4) 


■/iK^ V 1 / Ji(aft)Ji(afe) 
Jo{a(,i)) \ Jo(a^i) Jo(«^3) 


x sin(i?i + $3 — $4) > e* 1 ® 4 


Ji(q£s) 

Jo(a£s) 


Ji(o , $5) 

Jo(o^5) 


^ Jl (a^2) Jl(pC^ 3 ) 
yJo{a&) Jo (0^3) 
Ji (a^i) Ji(q^4) ^ 
Jo(oCi) Jo(a?4) y 


COS ($2 + $3 — $ 5 ) 


cos($i + $4 — $5) 


j / Jl( a C5) 

V Jo(a^5) 


4 ^*})-**+*-« 

^ Ji(q^i)Ji(q^4) \ 

v Jo(a£i) Jo(a^4) J 


x sin($i + i?4 — $5) > e ±lSs . 



